Preparing a highly degenerate Fermi gas in an optical lattice 
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We propose a method to prepare a sample of fermionic atoms in a three-dimensional (3D) optical 
lattice at unprecedentedly low temperatures and uniform filling factors. The process involves adia- 
batic loading of atoms into multiple energy bands of an optical lattice followed by a filtering stage 
whereby atoms from all but the ground band are removed. Of critical importance is the use of a non- 
harmonic trapping potential, taken here to be the radial profile of a high-order Laguerre-Gaussian 
laser beam, to provide external confinement for the atoms. For realistic experimental parameters, 
this procedure should produce samples with temperatures ~ 10 -3 of the Fermi temperature. This 
would allow the investigation of the low-temperature phase diagram of the Fermi-Hubbard model 
as well as the initialization of a high-fidelity quantum register. 

PACS numbers: 03.75.Ss, 32.80.Pj, 03.67.Lx, 37.10.De, 37.10.Jk, 05.30.-d 



Investigations of degenerate Fermi gasses loaded into 
optical lattices have indicated that these systems are 
ideal for creating a robust quantum register for quan- 
tum computing applications [l|, 12J as well as providing 
a testing ground for paradigm models of condensed mat- 
ter physics. Models currently under investigation include 
studying Fermi surfaces and band insulator states [3J, 
fermionic superfluidity in a lattice |4|, and transport 
properties of interacting fermions in one and three dimen- 
sional optical lattices [5|, [|| . These seminal experiments 
demonstrate the high precision and versatility available 
in simulating solid state systems with fermions in optical 
lattices. 

Theoretical studies of such systems have predicted that 
a number of exotic phases emerge at low temperatures, 
including quantum magnetic ordering and possibly d- 
wave superfluidity [7|, |8|]. However, temperatures low 
enough to observe exotic phases such as these are dif- 
ficult to achieve when optical lattices are loaded with a 
harmonic external confining potential. It has been the- 
oretically predicted [9j and experimentally observed [3j 
that fermions adiabatically loaded into an optical lattice 
with harmonic external confinement experience heating 
for all but very high initial temperatures and filling fac- 
tors (number of atoms per lattice site) [lOj] . 

Alternative methods to prepare fermionic atoms in op- 
tical lattices at low temperatures and/or high uniform 
filling factors include: cooling by adiabatic loading into 
a three-dimensional (3D) homogeneous trapping poten- 
tial with high filling factor [llj , defect filtering in a state 
dependent optical lattice [l| , adiabatic loading [2j and fil- 
tering [I2j of high entropy atoms from a ID lattice with 
harmonic confinement. 

In this Letter, we propose a method to prepare a highly 
degenerate Fermi gas in a 3D optical lattice using a box- 
like potential for external confinement and taking advan- 
tage of the Pauli exclusion principle to selectively remove 
atoms from multiply-occupied lattice sites. Specifically, 
we assume that the radial profile of a blue-detuned, high- 



FIG. 1: (color online) We filter high entropy atoms from a 
combined box-like potential (blue dashed line) and optical 
lattice (red solid line) by selectively removing atoms from 
all but the ground energy band. Amplitude modulation of 
the lattice potential can selectively transfer these atoms to 
high-lying bands via a two-photon transition where they then 
tunnel out of the region. Dramatic cooling results when the 
Fermi energy (prior to filtering) lies within the second band. 



order Laguerre-Gaussian (LG) laser beam provides con- 
finement along each cartesian axis. The atoms are pre- 
pared via a two step process: (I) adiabatically loading 
atoms initially confined in the LG trap into a superim- 
posed optical lattice, followed by (2) irreversibly filtering 
atoms from all but the ground energy band (see Fig. [I}. 
We find that when the Fermi energy of the system is 
sufficiently large, such that atoms begin to significantly 
populate the second energy band prior to filtering, con- 
siderable cooling is achieved; whereas for lower filling fac- 
tors heating is observed. 

The energy spectrum of a system of ultracold atoms 
is greatly affected by the addition of a 3D cubic optical 
lattice which can be formed by three perpendicular sets 
of retroflected Gaussian laser beams detuned far from 
resonance. In a homogeneous trapping potential, the lat- 
tice breaks the translational symmetry of the system, re- 
sulting in a series of discrete energy bands whose width 
and energy spacings are dependent on the intensity and 
wavelength of the laser light (see Fig. HJa)). Bezett and 
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FIG. 2: (color online) (a) Band structure for a 5 Er deep 
lattice, (b) Atoms can be selectively excited from the n to 
n + 2 energy band by modulating the amplitude of the op- 
tical lattice light. The excitation energy between bands 1 
and 3 (blue) and 2 and 4 (red) for fermions in all crystal 
momenta states are shown. As the depth of the optical lat- 
tice approaches 35 En, the excitation energy bands separate, 
thereby demonstrating the feasibility of performing the band 
selective excitations required in our filtering method. 



Blakie demonstrated that this energy band structure can 
be exploited to dramatically increase the degeneracy of 
the sample for a homogeneous system 



111. For a dense 



atomic gas, with a filling factor greater than unity, ap- 
plication of the lattice increases the Fermi energy, since 
it lies within the first excited band, and compresses the 
Fermi surface resulting in a dramatic reduction in the 
degeneracy temperature T/Tp, where Tp is the Fermi 
temperature. 

This band structure also permits state-selective opera- 
tions to manipulate and probe the energy distribution of 
the sample. One such method involves modulating the 
depth of the optical lattice to selectively excite atoms 
from the n to n + 2 energy band with no change in the 
crystal momentum q [13J. In contrast to non- interacting 
Bose systems, where a single q can be macroscopically oc- 
cupied, a Fermi system necessarily begins to fill the bot- 
tom band and Aq — transitions must be excited for all 
occupied values of q. In Fig.^b), we show band excita- 
tion energies as a function of lattice depth for n = 1 — v 3 
and n = 2 — » 4 transitions spanning all q within a Bril- 
louin zone. By loading the sample into an optical lattice 
with a depth of Vq — 35_Er (where En = K k 2 /2m is 
the recoil energy and fc is the wavenumber of the laser 
light), we find that these transitions are well resolved. 
It is therefore possible to apply a filtering process which 
selectively removes atoms from all but the ground energy 
band. Using adiabatic rapid passage, population may be 
selectively transferred from n — 2 — > 4 by sweeping the 
amplitude modulation frequency from below to above all 
2^4 transition frequencies while remaining below the 
lowest 1 — > 3 transition frequency. Then, lowering the 
height of the trapping potential allows atoms in the third 
and higher energy bands to tunnel out of the system. 

In order to experimentally approximate the homoge- 
nous lattice potential described above, we consider the 



addition of a box-like external potential produced by a 
blue-detuned, t -order, LG laser beam with a radial pro- 
file 
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at the beam waist wq. For a given charge £, the peak 
value Vpeak of the potential occurs at r max = wo\fl/2 
and the width of this peak decreases with decreasing wq. 
Therefore, for a given trap size r max , the LG profile more 
closely approximates a box potential when u> is reduced 
and I is correspondingly increased. Trapping of ultracold 
gases has been demonstrated in single or crossed beam 
configurations of LG beams up to I = 16 [14J. Il5lll6l|. 

Along a given cartesian axis, we take the single particle 
Hamiltonian to be 



H(x) 
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2m dx 2 
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+Vq cos (kx + 4> x ) + -^fnuj x , 



where the third term represents a lattice potential of 
depth Vq and phase offset <f> x . We also include a har- 
monic term that arises if red-detuned Gaussian beams 
are used to produce the lattice potential; in this case 
oj oc \/Vo- The ID eigenvalues and eigenfunctions for a 
given depth of the optical lattice are calculated by nu- 
merically diagonalizing the Hamiltonian (Eq. ^ using 



the method described in [17|. For sufficiently shallow 



lattice depths, the low energy eigenstates are delocalized 
and closely approximate Bloch states in the first band of 
a homogeneous system. However, higher energy states 
are either localized at the edges of the trap (i.e. near 
x = r max ) or delocalized and correspond to Bloch states 
in higher bands. While a band structure picture is not 
strictly valid for this inhomogeneous system, we classify 
the set of eigenfunctions without nodes to constitute the 
first band, ei ba nd,n- 

We extend this model to three dimensions by assuming 
a separable Hamiltonian 7?3D = H(x) + H(y) + H(z). 
For simplicity, we assume equal lattice depths in each 
direction. The 3D spectrum (£ m ) for a given depth of the 
optical lattice is then generated by calculating all possible 
combinations of the sum £ m = e-i + tj + £& for all values of 
the ID eigenenergies (e p ) in each spatial direction. The 
3D energy spectrum for energy states in the first band 
of the optical lattice is calculated in a similar manner 

Riband, m — ^lband,i ~r Clband.j t ^lband,fcj- 

In calculating thermodynamic quantities during the 
proposed cooling method, we assume constant thermal 
equilibrium before and after the selective removal of 
atoms from high- lying bands. Equilibrium is main- 
tained by elastic collisions in a 50/50 mixture of spin-1/2 
fermions and changes in the trapping potential are adi- 
abatic with respect to the rethermalization time scale. 



However, we also assume that the interactions are weak 
enough to not significantly modify the single-particle en- 
ergy spectrum, £ m . We therefore use £ m when calculat- 
ing the following quantities: 

^ 1 
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exp[(£ m - jj)/k B T] 
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exp[(£ m - y)/k B T] 



— = 2^1n[l + exp[(£ m -/!)/fc B T]] 
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where T is the temperature, /x is the chemical potential 
of an atom in either spin state, N is the total number of 
atoms, E is the total energy in the system, and 5* is the 
total entropy. The degeneracy temperature is given by 
T/T F where k B T F = £ N , the energy of the N th eigen- 
state of the multi-band, 3D spectrum. After filtering, 
the thermodynamic quantities in Eq. [3] are calculated for 
atoms only in the first band using £iband,m- 

Our proposed method for cooling the atoms is com- 
prised of (1) an adiabatic increase in the lattice depth 
starting from zero, (2) a non-adiabatic selective filtering 
of atoms and (3) an optional adiabatic change to a final 
lattice depth. To calculate changes in thermodynamic 
quantities during these stages we use the following pro- 
cedures. For adiabatic changes of the potential we (1) 
calculate S for a given N and initial temperature 2] us- 
ing the energy spectrum for the initial potential and (2) 
numerically solve for /if and the final temperature Tf in 
Eq. [3] using the spectrum for the final potential, assum- 
ing N and S are conserved. In contrast, for the selective 
filtering stage we (1) start from a thermalized sample of 
N[ atoms at temperature 7] for a given spectrum, (2) cal- 
culate, given this distribution, the energy E[ and number 
N{ for atoms restricted to the first band, and (3) solve for 
Hi and the temperature Tf in Eq. [3] using the multi-band 
energy spectrum assuming the sample equilibrates with 
total energy E{ and number JVf. 

We consider a 50/50 spin mixture of 6 Li atoms ini- 
tially trapped in a LG trapping potential with £ = 12, 
V pca k — 35 Er and r max = 13.5 ^m. For reasonable 
lattice beam properties (k = 27r/1064nm and a waist 
of 200 (im) we find uj = 2ir (586Hz) for the final lattice 
depth Vo,f = 35 _Er. The final degeneracy temperatures 
after adiabatic loading and filtering, along with the fi- 
nal atom number are shown in Fig. [3] for various initial 
degeneracy temperatures and sample sizes. In each case 
4>x = 4>y — 4>z = 0- This data shows that the thermo- 
dynamic properties of the system are highly dependent 
on the initial filling factor and can be separated into two 
distinct regions A and B. The vertical dashed line which 
separates the regions represents the number of atoms at 
which the Fermi energy enters the second band. 

In region A, the Fermi energy before filtering lies be- 
low the second energy band. For very low filling factors, 
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FIG. 3: (color online) As a function of initial atom number 
we report (a) the final degeneracy temperature and (b) the 
final atom number after implementing the proposed cooling 
and filtering procedure for various initial temperatures be- 
tween 0.005 and 0.3 Tf- The vertical dashed line represents 
the number of atoms for which the Fermi energy enters the 
second band. The trap and lattice parameters are as described 
in the text. 



an increase in T /Tp is observed for all initial tempera- 
tures. This increase in T/Tp occurs because Tp decreases 
more than T as the lattice depth increases. Additionally, 
Fig. E3b) demonstrates that we are not significantly fil- 
tering atoms for low initial temperatures. As the Fermi 
energy approaches the second energy band, we see a dra- 
matic decrease in the final T/Tp. In this regime, where 
atoms are beginning to significantly occupy the localized 
energy states at the edge of the trapping potential, a 
dramatic increase of the Fermi energy is observed. 

In region B, the density is such that the Fermi energy 
before filtering lies within the second band. In this re- 
gion the adiabatic increase of the lattice depth results in 
a dramatic increase in Tp, a substantial reduction in the 
temperature T 11(, and allows for a significant reduc- 



tion in entropy during the filtering stage. We find that 
significant cooling is achieved for initial temperatures in 
the vicinity of T — 0.1 Tp. Above this initial tempera- 
ture cooling is less efficient. Below this initial tempera- 



ture, the final T/Tp after filtering saturates. As can be 
seen in Fig. [3jb), N{ is extremely insensitive to fluctua- 
tions in N[ for low initial temperatures. For example, at 
3] = 0.05 T F , a variation of ±10% around JV ; = 1.6 x 10 5 
yields a variation of only +0.09%/ - 0.2% in N f . 

The cooling efficiency and number filtering were de- 
pendent on the choice of phases 4> x , 4> y and 4> z due to 
the sensitive effect these phases had on the the location 
of localized edge state eigenenergies relative to the Fermi 
energy. To study this effect, we modeled the system al- 
lowing the phase in each direction to be independently 
selected from the set <p a = (0, ir/10, ...,ir/2). We consid- 
ered samples with an initial temperature 7] = 0.05 Tp 
and an initial number iVj = 1.6 x 10 5 atoms, parameters 
within the saturated regime for all choices of phase and 
close to optimal for cooling. From the set of all possible 
phase combinations, we find an average final tempera- 
ture 7f = 0.0031 Tp where 10% of the ensemble achieved 
temperatures below T w = 0.0023 Tp and 90% were below 
Tgo = 0.004 Tp. From this same set, we find an average 
final number N t = 1.20 x 10 5 , with N 10 = 1.18 x 10 5 , 
and iVgo = 1.22 x 10 5 . The filtering process further re- 
sults in a substantial reduction in entropy. The initial 
entropy per atom Si = 0.28 fee is reduced to an average 
final value of S{ — 0.024 fc^, with Sio = 0.014 ks, and 
sgo = 0.033 fcs. 

It is in general possible to prepare atoms at a low 
T/Tp in a shallow lattice potential, if so desired, by adi- 
abatically reducing the lattice depth after the filtering 
stage. Continuing the example from above, when the 
lattice depth is reduced to 5 En we find an average fi- 
nal temperature T f = 0.002 T F , T w = 0.0013 T F , and 
T 9Q = 0.0028 T F . 

We now consider the effects of the charge £ of the LG 
beam for samples with an initial 2] = 0.05 Tp, phase 
(j) x = (j) y = 4> z — 0, final lattice depth of 35 En, and var- 
ious initial atom numbers. For each ^-value, the waist 
of the LG beam is adjusted such that the number of 
states below the second energy band is held constant at 
1.22 x 10 5 . As shown in Fig. [4] the cooling efficiency of 
this procedure is highly dependent on the charge. Note 
that for t = 1, which approximates harmonic exter- 
nal confinement, the final degeneracy temperature T{/T F 
never drops below its initial value of T-JT F = 0.05. For 
£ > 8, the minimum degeneracy temperature saturates 
to Tf/Tp < 0.003. For higher values of £, the extent of 
the saturation regime grows. We believe that this satu- 
ration is caused by localized atoms at the edges of the 
LG potential rethermalizing into higher energy bands. 

In this Letter we proposed a method for preparing a 
sample of highly degenerate fermions by adiabatic load- 
ing into a combined optical lattice and "box-like" trap- 
ping potential followed by selective removal of atoms 
from all but the ground energy band. Numerical calcu- 
lations for sample sizes ~ 10 5 predict that temperatures 
~ 10~ 3 Tp can be prepared in this manner. This method 
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FIG. 4: (color online) The effects of the charge £ of the 
Laguerre - Gaussian trapping potential on the efficiency of 
our proposed cooling and filtering method. For each data 
set, the initial temperature 7] = 0.05 Tf and the phases 
4>x = 4> y — 4>z = 0. For all £ values, the number of atoms 
at which the Fermi energy enters the second band (vertical 
dashed line) is held constant. 



is robust against initial number and temperature fluctu- 
ations for a sufficiently cold initial sample of atoms and 
yields samples with little variance in the final number. 
While the selective removal of atoms must occur in a 
deep lattice (in order to spectrally resolve the band exci- 
tations), subsequent reduction of the lattice depth, if de- 
sired, yields a modest amount of additional cooling. We 
expect that this method can be scaled to larger samples 
for which still lower degeneracy temperatures would be 
attained due to the diminished role localized edge states 
would play. Further, the "box-like" trapping potential 
offers an ideal spatial profile for simulating solid state 
physics with degenerate atoms in optical lattices as the 
relatively flat central region allows for a large number 
of delocalized states while the curvature at the edges of 
the traps removes the constraint of loading exact atom 
numbers to realize insulating states. Atoms prepared in 
this manner should be sufficiently cold to explore quan- 
tum spin phases of fermionic atoms which are currently 
inaccessible, and could provide a physical realization of 
an essentially perfect quantum register. 
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